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We present the conditional time evolution of the elec- 
tromagnetic field produced by a cavity QED system in the 
strongly coupled regime. We obtain the conditional evolution 
through a wave-particle correlation function that measures 
the time evolution of the field after the detection of a pho- 
ton. A connection exists between this correlation function 
and the spectrum of squeezing which permits the study of 
squeezed states in the time domain. We calculate the spec- 
trum of squeezing from the master equation for the reduced 
density matrix using both the quantum regression theorem 
and quantum trajectories. Our calculations not only show 
that spontaneous emission degrades the squeezing signal, but 
they also point to the dynamical processes that cause this 
degradation. 

OCIS codes: (270.0270) Quantum Optics; (270.6570) 
Squeezed States; (270.2500) Fluctuations, relaxations, 
and noise. 



I. INTRODUCTION 

Squeezing experiments focus on the control and study 
of quantum fluctuations. In a squeezed state there is a 
redistribution of the fluctuations between the different 
quadratures of the electromagnetic field. The modifica- 
tion is measurable through the spectrum of squeezing 
and its comparison with that of the vacuum state. In 
this paper we analyze the time evolution of the quantum 
fluctuations of the electromagnetic field, the complemen- 
tary picture of the spectral density studied in squeezing. 
We use a new formalism based on a third order correla- 
tion function of the field. This wave-particle correlation 
combines conditional measurements with the process of 
quantum measurement. The quantum fluctuations are 
related to the measurement process and, in a strongly 
coupled system, those fluctuations can be orders of mag- 
nitude larger than the steady state. 

By combining the wave and particle aspects of light, 
this correlation function brings a new theoretical and 
experimental [^|j3j understanding to the study of the 
quantum properties of squeezed light. This is because 
the wave-particle correlation function and the spectrum 
of squeezing form a Fourier transform pair when third 



order moments of the field fluctuations are neglected. 
We can now analyze the conditional time evolution of 
the fluctuations that give rise to the redistribution of the 
noise between the two quadratures of the electromagnetic 
field. 

Work on the spectrum of squeezing is a well estab- 
lished area of research The interaction between a 
single cavity mode of the electromagnetic field and sev- 
eral two level atoms produce non-classical fields which 
exhibit squeezing. Work on this system started with the 
problem of optical bistabihty (OB) ||-|lO|. These OB 
systems are treated in the limit of small noise, such that 
the fluctuations are only a minimal perturbation on the 
steady state. More recent work has focused on a regime 
where there is no small size parameter and fluctuations 
dominate the behavior of the system This general 
area is cavity QED [|l2|. While there are experimental 
realizations in the microwave and optical regimes, our 
discussion will only consider the latter. In both systems 
the reversible coupling between the atoms and the cavity 
mode is larger than the irreversible loss of coherence from 
spontaneous emission or cavity decay. A number of op- 
tical studies, experimental and theoretical, have focused 
on the non-classical features as observed in the intensity 
correlations, eg. photon antibunching [p"3|-p^. Our recent 
observations look, however, at the time dependence of the 
conditional fluctuations of the field |^ and its connection 
with the spectrum of squeezing. 

We found from experiment that the spectrum of 
squeezing was degraded when the conditional fluctuation 
oscillated around a value different from the steady state 
field amplitude. In order to understand this behaviour we 
apply the time domain formalism to the strong coupling 
regime of cavity QED. To keep the discussion manage- 
able we only treat the case of one or two atoms fixed 
in space. The work follows the general philosophy es- 
tablished in Refs. where cavity QED is analyzed 
without a small parameter. 

The paper is organized as follows. In section II, we 
describe the wave-particle correlation function and its 
connection to the spectrum of squeezing. Section III dis- 
cusses the properties of the cavity QED system. Section 
IV shows the different techniques utilized to calculate the 
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spectrum of squeezing and the wave-particle correlation 
function. We present our results in section V, the discus- 
sion of those results in section VI, and our conclusions in 
section VII. 



II. THE WAVE-PARTICLE CORRELATION 
FUNCTION 

Figure 1 shows the apparatus that measures the cor- 
relation between the amplitude and the intensity of the 
electromagnetic field. The particle (photon) produces a 
trigger 'click' in an avalanche photo diode (APD) and 
conditionally the wave (electromagnetic field amplitude) 
gets recorded in the photocurrent output of a balanced 
homodyne detector (BHD). The observations of Ref. ||] 
consider a conditioned field evolution described by a 
third order correlation function of an electromagnetic 
field mode b (signal) that has a non-zero, steady state 
average (b) = A. The measured correlation is 
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where Qe = {bexp{—i9) + b''exp{i6))/2 is the quadrature 
amplitude selected by the local oscillator phase 9 and rj 
is the coupling efficiency into the BHD. Carmichael et 
al. have generalized Eq. (||) to any source , but for the 
purpose of this paper it is sufficient to limit the discussion 
to the detection scheme of Fig. 1. 

We separate the fluctuations from the mean field writ- 
ing. 



b = Aexp(i6') + Ab 



(2) 



The correlation function can be rewritten in terms of the 
field fluctuations. In the limit of Gaussian fluctuations we 
may neglect third-order moments of the field fluctuations 
and the correlation function is iBl 
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where AQo" = (A6 + A6^)/2 and ^(r) is the residual shot 
noise with the following correlation function 



-exp(— Ft) 



(4) 



where F is the BHD bandwidth. It is apparent from 
Eq. (^ that once the conditioned field is normalized for 
a fixed then the shot noise, ^{t), is determined 

by the detection bandwidth and the number of samples 
averaged. Therefore one must average over many start 
"clicks" to obtain a good signal to noise ratio. 
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FIG. 1. Schematic of the wave particle correlator. 

The spectrum of squeezing is the cosine transform of 
the normalized correlation function, 



S{d, i^)^4F dr cos(27ri/r) \h^i 
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where F — 2n{b%) is the photon flux into the correlator 
and hg{T) is the average of /ie(r) with respect to ^(t). 
This average is independent of the fraction of light, 1 — r, 
sent to the BHD. 

The wave-particle correlator has two practical advan- 
tages over a standard squeezing measurement. First, the 
fluctuations in time reveal the entire spectrum of squeez- 
ing up to a frequency set by the detector bandwidth. The 
second is that the BHD efficiency only enters through 
the residual shot noise in Eq. (1). This shot noise can 
be made to vanish by averaging over many trigger starts. 
This implies that the spectrum of squeezing is not de- 
graded by imperfect detector efficiencies and propagation 
errors which beset the standard squeezing measurement. 



III. CAVITY QED SYSTEM 

The fiuctuations in the signal beam in Fig. ^ must be 
squeezed by a non-linear source. The source we consider 
here originates from a cavity QED system. The cavity 
QED system consists of a single mode of the electromag- 
netic field interacting with a collection of two-level atoms 
(see Fig. 2). Two spherical mirrors that form a standing- 
wave optical cavity define the field mode. We consider 
only a single or a pair of two level atoms optimally cou- 
pled to the mode of the cavity. Dissipation plays an im- 
portant role as both the atoms and the field couple to 
reservoirs. An atom can spontaneously emit light into 
modes other than the preferred mode, and light inside 
the cavity can escape through the mirrors. We assume 
that the fractional solid angle subtended by the cavity 
mode is small enough so that we do not have to make 
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corrections to the atomic decay rates. A coherent field, 
injected through one of the mirrors, drives the system. 
The signal beam in Fig. [l] is the light that escapes from 
the cavity through the output mirror. 
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FIG. 2. Schematic of the cavity QED system. 

The Jaynes-Cummings Hamiltonian describes the 
interaction of a two-level atom with a single mode of the 
quantized electromagnetic field. 
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where a± and are the Pauli spin operators for rais- 
ing, lowering, and inversion of the atom, and , a are the 
raising and lowering operators for the field. The eigen- 
states for Eq. (||) reveal the entanglement between the 
atom and the field. The spectrum has a first excited 
state doublet with states shifted by ±g from the uncou- 
pled resonance. The dipole coupling constant, g, is given 
by: 
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where fj, is the transition dipole moment, lu is the tran- 
sition frequency, and V is the cavity mode volume. The 
field decays out of the cavity at the rate of k and, consid- 
ering only radiative decay, the atomic inversion decays 
at the rate of 7 = 1/t (t is the radiative lifetime of the 
atomic transition). 

From g, 7 and k we can construct two dimensionless 
numbers from the OB literature that are useful for char- 
acterizing cavity QED systems: the saturation photon 
number uq and the single atom cooperativity Ci. They 
scale the influence of a photon and the influence of an 
atom in the system. These two numbers relate the re- 
versible dipole coupling between a single atom and the 
cavity mode (g) with the irreversible coupling to the 
reservoirs through cavity (k) and inversion decays (7) 
by 



(8) 



and for a plane wave cavity field with the atoms sitting 
at the anti-nodes, 



no = 7'/%' 



(9) 



They define the different regimes of OB and cavity 
QED. The strong couphng regime of cavity QED, no < 1 
and Ci > 1, implies large effects from the presence of a 
single photon and of a single atom in the system. Con- 
sequently, a change of one photon is a large fluctuation 
in the system. 

We only consider the resonant case here, with the driv- 
ing field on resonance with the atomic transition and cav- 
ity mode frequencies. This means that we are only able 
to explore dynamical changes in one quadrature. In or- 
der to generalize to other phases one must either drive 
the system off resonance or add a coherent offset like the 
one discussed in Ref. 

Two dimensionless fields and intensities follow from 
the OB literature that allow us to make contact with 
experiments: The intracavity field (intensity) with atoms 
in the cavity x = {a)/^/no; {X = (a^a)/no), and the field 
(intensity) without atoms in the cavity y = £ j k^Jtiq^ 
{Y = j/^) where 2£^/k is the input photon flux. 



IV. CALCULATIONS OF THE WAVE-PARTICLE 
CORRELATION FUNCTION IN CAVITY QED 

We present the two different methods used to calculate 
the wave-particle correlation function. The first is a di- 
rect calculation from the master equation and the second 
is a quantum trajectory simulation. 



A. Direct Calculation from the Master Equation 

A master equation of the Lindblad form can be derived 
for the system shown in Fig. ^ using standard techniques 
(see for example j2^]). This equation for the reduced 
density operator p for the case of A^-atoms identically 
coupled to a single cavity mode reads 

p = Lp — E\a^ — j}\ ^ 5[a^S'_ — a<S'_|_, p\ (10) 
+ K,(2apa' — a} hp — pa) a) 

N 

where S± — (Tj_ are the collective raising and lowering 
operators for the atoms. 

Early attempts to solve this master equation are found 
in the OB literature §,0]. They involved converting the 
master equation into a c-number Fokker-Planck equation 
via a large N expansion. A connection was then made to 
a set of Ito stochastic differential equations. This master 
equation could thus be solved analytically, but only for 
the case of small fluctuations. 

Our cavity QED system operates in the strongly cou- 
pled regime where the effective number of atoms and the 
size of the quantum fluctuations, relative to the steady 
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state, will not permit us to make any of the assumptions 
mentioned above. Instead we numerically solve the mas- 
ter equation in a manner similar to that described by 
Brecha et al . We transform the master equation into 
frequency space and apply the quantum regression the- 
orem of Lax 1^ to calculate the spectrum of squeezing 
directly. We check the accuracy of our results by increas- 
ing the maximum number of photons included in the ba- 
sis, rimax, until the value of the squeezing at any given 
frequency is accurate to three significant digits. This ap- 
proach permits us to calculate the spectrum of squeezing 
for the cavity QED system away from the low intensity 
limit and it allows us to consider the case of either one 
or two atoms. 



B. Quantum Trajectories 

Studying quantum trajectories with the wave-particle 
correlator requires some justification. This section pro- 
vides that justification. We begin with a review of some 
of the fundamental concepts behind quantum trajectory 
theory. Much of the important work done in the field of 
quantum trajectories and stochastic Schrodinger equa- 
tions can be found in the review article by Plenio and 
Knight and references therein. This presentation fol- 
lows the work of Carmichael [po|| . 

An alternative approach to studying a master equation 
of the Lindblad form can be formulated with quantum 
trajectory theory. To create quantum trajectories one 
must unravel the master equation. We are only interested 
in those unravellings which correspond directly to an ex- 
perimental measurement. This measurement probes the 
system by opening channels between the system and the 
environment. Excitations into these channels are repre- 
sented by the following superoperator 



(11) 



where fc is a function of the creation and annihilation 
operators, sj and Si, for the i*'* system operator coupled 
to the environment. The probability for the system to 
decay through one of these channels depends on the oc- 
cupation of that particular mode, 



(12) 



Quantum trajectories also provide a picture of what is 
happening to the system between measurements. Equa- 
tion (11) describes how measurements are performed on 
the system. To see what happens between measurements, 
subtract Eq. (0) from Eq. (|l^). 



(13) 



If we assume that the system state is initially pure then 
Eqs. ( pT| ) and (|l^) guarantee that the conditional density 



operator can be written in the factorized form p = \'ip){ip\. 
Therefore Eq. ( |l3| ) implies that a propagator exists for 
the system wavefunction between measurements. This 
propagator is a modified, non-unitary, Hamiltonian with 
a real and imaginary part, -ffsys = Hn + iHi. 

Below is the quantum trajectory algorithm which prop- 
agates the system wavefunction forward in time. 

1. Choose an initial state for the system. 

2. Calculate the probability for the system to decay 
through each of the loss channels from Eq. (|l^) . 



P, = {k\si)k{h))dt 



(14) 



3. Associate a uniformly distributed random number 
between and 1 with each loss channel in step 2. If the 
probability for the channel is greater than its corre- 
sponding random number then collapse the system wave- 
function. 



m + dt)) = k{h)m)) 



(15) 



In the unlikely event of multiple collapses, choose only 
one collapse via another random number. 

4. If the probability for loss through each channel is 
less than the corresponding random numbers, then prop- 
agate the system wavefunction with the effective Hamil- 
tonian from Eq. (13), 



m+dt)) = {i 



ifi 



''-dt)\m) 



(16) 



5. Normalize the system wavefunction and repeat from 
step 2. 

This is the basis for quantum trajectories. The master 
equation describes the evolution of the reduced density 
matrix. A connection is made between this density ma- 
trix and a system wavefunction. This wavefunction leads 
to measurements via losses through the system's decay 
channels. These measurements update our knowledge of 
the system wavefunction. This allows one to study how 
measurements performed on a system will, through con- 
ditioning, affect the system state. 

In order to unravel Eq. ( [l0| ) with the wave-particle 
correlator we need to first add the BHD to the system. 
The BHD consists of a local oscillator which is modeled 
by a driven, single mode coherent field. We write the 
master equation for the system plus local oscillator mode 
as 

p = £[a^ — a, p] + g[a) S- — aS+, p] (17) 
+ K{2apa) — a) ap — pa) a) 
+ KLoi'2c0 ~ c^cp - 0c) + klox[c^ - c, p] 

N 

+ -il2^(2a'_pa\ - cj\alp - pa^al), 

where & and c represent the raising and lowering oper- 
ators for the local oscillator mode. The coherent field 
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occupies a mode whose strength we denote by x- Under 
a BHD scheme the beam sphtter combines the signal and 
local oscillator fields equally {R = 1/2), as in Fig. 1, to 
give the following total measured fields at photodetectors 
1 and 2 



BHD 



,2 = iz-^KLoc + \/ k(1 — r)a 



(18) 



while the measured field at the photon counting detector 
gives 



(19) 



We assume that the signal mode does not affect the lo- 
cal oscillator mode. This implies that the density matrix 
separates into a piece which corresponds to the cavity 
QED system and a piece which corresponds to the local 
oscillator, 



Ps\x) (xl- 



(20) 



We define the local oscillator flux in terms of the 
strength of the local oscillator mode and its cavity de- 
cay rate. 



f = KLolxl' 



(21) 



From Eq. ( |1 1| ) we construct the master equation terms 
which correspond to Eqs. ( p^p^ ) and spontaneous emis- 
sions. 



5bhdi.2/5. = (iV/e'*' + v/2k(1 - r)a) 

xps{±^/fe-''' + ^2Kil-r)a^) (22) 



Scount Ps = 2Krapsa 



5, 



pontPs 



(23) 



(24) 



After subtracting Eqs. ( p2[p4| ) from Eq. ( |17|), we arrive 
at an expression which corresponds to Eq. (|l3|). 

(^ — '5bhD1-2 ~ ^count — Sspont )Ps = f [o^ — &, Ps] 

+g[a^S- — dS+, ps] — K{a^aps + Pso'a) 



N 

-l/'^^iK^-Ps + Ps^i^-) - fPs 



(25) 



Equation ( p5| ) provides a modified Hamiltonian which 
propagates the system's wavefunction between measure- 
ments. Equations ( p^ ) and (^3|) govern how "clicks" at 
the BHD and APD disrupt that evolution. 

There are two types of detections, not including spon- 
taneous emissions, which occur under the conditional 
BHD measurement. The first comes from Eq. ( p^ ) 
which corresponds to detections at the APD. According 
to Eq. (12) these collapses occur with a probability of 



2nr{a) a)dt and they produce a signifigant change in the 
wavefunction as is evident from Eq. (|l5|). 

The second type of collapse corresponds to the "clicks" 
at the BHD. The strength of the local oscillator tells us 
that there are many BHD "clicks" within a time interval 
set by K^^. Eq. ( p5| ) shows that the effect of each of these 
"clicks" on the system wave function is almost negligible. 
This implies that a straight forward trajectory simulation 
with the algorithm described above in Eqs. ( p^[]l6| ) will 
not provide an efficient method for modeling the wave- 
particle correlator. Instead we may associate a photocur- 
rent with the detections registered by the BHD. The su- 
perposition of the cavity field with the local oscillator 
field in Eq. ( |l8| ) implies that the BHD photocurrent and 
the system wavefunction evolution are not independent. 
With the principles outlined in Sects. 8.4, 9.2, and 9.4 of 
pO| we combine Eqs. ([2^p5|) to simultaneously calculate 
this photocurrent from the BHD setup and then evolve a 
corresponding stochastic Schrodinger equation forward in 
time. This analysis leads to the following expressions for 
the difference in photocurrents between photodetectors 1 
and 2 in Fig. 1, and wavefunction propagation between 
photodetections at the APD and spontaneous emissions 
out the sides of the cavity. 



di 



-T{idt - ^j8K{l-r){de)cdt + dWt) 



(26) 



dm^ = \-^dt + ^2k(1 - r)aexp(-i6l) 
in 



- r){de)cdt + dWt)\^)c (27) 

where F is the BHD bandwidth and dWt is the same 
Wiener noise increment in both equations. 

This BHD difference current, z(t), along with a set of 
start times {tj}, form a stochastic measurement record. 
The source quasimode a is in a quantum state IV'rec (O) 
conditioned on the past record. We simulate Eqs. ( Pq - 
p7| ) on a computer. By sampling an ongoing realization 
of i{t) for many "start" times (APD detections realized 
concurrently) we may calculate the following averaged 
photocurrent. 



H(r) 



(28) 



A connection exists between the averaged photocur- 
rent and the wave-particle correlation function ||l[. In 
the limit of large bandwidth this connection reads 



heir) 



n{T) 



{d)^8K{l-ry 



(29) 



The dual nature of the measurement process provides 
one of the many strengths of the conditional field mea- 
surement. It allows the use of quantum trajectory theory 
to unravel the master equation in two distinct ways. A 
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numerical simulation of Eqs. (]26|-p7D reproduces the ex- 
perimental results for the averaged photocurrent. We can 
understand some of the mechanisms which lead to this 
averaged result by simulating the algorithm described in 
Eqs. and replacing Eq. ^ with Eq. @. Set- 

ting r sa 1 reduces the conditional homodyne simulation 
to a cleaner photocounting simulation. This leads to an 
understanding of the dynamical processes which create 
the spectrum of squeezing. 



V. RESULTS 

The simplest parameter to change during the course of 
the experiment is the strength of the driving field. This 
motivates our study of the effect of driving intensity on 
the spectrum of squeezing. The strength of the driving 
field is parameterized in terms of the saturation photon 
number, tIq. Therefore weak fields imply X <^ 1, for 
intermediate fields X « 1, and strong fields correspond 
to X ^ 1. We begin by reporting the results of a low 
intensity calculation with both the quantum trajectories 
and the direct calculation from the master equation. We 
then increase the intensity of the driving field and also 
the number of atoms in the cavity from one to two. This 
degrades the spectrum of squeezing in a way which we 
hope to understand through the study of quantum tra- 
jectories. 



A. Low Intensity Results for /i(0°,r) and 5(0°, i^) 

Figure ^ shows the equivalence of the third order cor- 
relation function ft,(0°,T) (i) and the spectrum of squeez- 
ing (ii) for a strongly coupled cavity QED system. The 
dashed line in Fig. |3|ii is the spectrum of squeezing cal- 
culated directly from the quantum regression theorem. 
The solid line is the Fourier transform (see Eq. (|^)) 
of Fig. H which comes from averaging the photocur- 
rent from a quantum trajectory simulation over 55000 
"starts" . Both approaches show the damped Rabi oscil- 
lations which follow and precede a photodetction. In the 
weak field excitation limit. Rice and Carmichael de- 
rived an analytical expression for the spectrum of squeez- 
ing which agrees with these results. 

Detecting a photon in the APD collapses the state of 
the cavity field. The oscillations in the field after a de- 
tection show the collapsed state oscillating back towards 
the steady state. In the Gaussian noise approximation 
the time symmetry of the measured correlation function, 
Eq. (0), guarantees the oscillation before the detection. 
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FIG. 3. i. Wave-particle correlation h{Q° ,t) for very low 
intensity excitation calculated from the quantum trajectory 
implementation of the conditioned homodyne photocurrent. 
The following parameters were used: (<;, ft, 7, r)/(2/7r) = 
(38.0, 8.7, 3.0, 100) MHz, X = 2.99 x 10"^, r = 0.5 and iVstarts 
= 55000. ii. Spectrum of squeezing calculated from the cosine 
Fourier transform of /i(0°, r). The dashed line is the spectrum 
of squeezing calculated directly from the master equation with 
the quantum regression theorem. 



B. From one atom to two atoms at intermediate 
intensity 

It has not been possible to obtain analytical expres- 
sions for the spectrum of squeezing in the case of interme- 
diate driving fields. In this regime one cannot linearize 
the system fluctuations about the steady state field as 
was done in the strong field case |^^. It is also impossi- 
ble to perform a small parameter expansion in the driv- 
ing field as was done in the weak field case ||l^ . We are 
instead forced to rely on numerical calculations to give 
us insight into a system's spectrum of squeezing. This 
approach has its limitations. 

The spectra of squeezing for increased excitation in 
the case of one and two atoms have distinct quantitative 
differences. Fig. ^ shows one such spectrum calculated 
for a driving field oi £ / k. — 1.50 and one atom. Fig. |^ 
shows the result for two atoms with £/k = 0.975. We 
calculated both plots with the same decay parameters 
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but different atom field couplings, g to keep the Rabi 
frequency, g\/N , the same. 
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FIG. 4. The spectrum of squeezing calculated directly from 
the master equation with the quantum regression theorem 
for N=l and high intensity. X = 104.0, (5, k, 7)/(27r) = 
(38.0,8.7,3.0) MHz. 
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FIG. 5. The spectrum of squeezing calculated directly from 
the master equation with the quantum regression theorem for 
N=2 and intermediate intensity. X = 1.36, (gr, k, 7)/(27r) = 
(38.0/^2,8.7,3.0) MHz. 

The width and shape of the central peaks are very dif- 
ferent. For the case of the single atom it is clear, from 
Fig. ^ and ||, that the cavity decay rate dominates the 
width. In order to see a peak at zero frequency one must 
drive the one atom system much harder than the two 
atom system, note that the intracavity intensity, normal- 
ized by the saturation photon number, is nearly 100 times 
greater for the one atom case. This can be understood 
from the point of view of an atom saturated with a strong 
driving field. At this point most of the photons that enter 
the cavity will bounce between the two mirrors without 
ever interacting with this single atom. This leads to less 
squeezing and a peak in the zero quadrature spectra at 
zero frequency with a width which is dominated by the 
cavity lifetime. 

When we plot the full-width at half-max (FWHM) for 



different driving fields in the two atom system. Fig. |^, 
we can see that spontaneous emission plays a more in- 
fluential role. This is not surprising because the cou- 
pling to the environment leads to a degradation of the 
squeezing signal, and these incoherent processes appear 
at zero frequency in the interaction picture. The ques- 
tion we hope to answer with quantum trajectories is the 
following: What is the mechanism by which spontaneous 
emission gives rise to this degradation in the squeezing 
signal? 




FIG. 6. Dependence of the width of the peak at zero fre- 
quency in the squeezing spectrum as a function of the in- 
put driving strength for N=l. The following parameters were 
used: {g,K,j)/2n = (38.0,8.7,3.0) MHz. The FWHM has 
been normalized by the cavity decay rate, k. 
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FIG. 7. Dependence of the width of the peak at zero 
frequency in the squeezing spectrum for three different 
spontaneous emission rates as a function of the input 
driving strength for N=2 and the following parameters: 
(5r,K,7i,72,73)/27r = (38/^2,8.7,3.0,1.0,0.5) MHz. The 
FWHM for each plot has been normalized by its respective 
spontaneous emission rate, 7i. The circles are for 7 — 0.5, the 
triangles are for 7 = 1.0, and the x's are for 7 = 3.0. Notice 
the difference in driving strengths, e/n, between this figure 
and the previous one. 
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C. High Intensity 

A large positive peak in the spectrum of squeezing ap- 
pears at zero frequency when we increase the strength of 
the driving field to the case when there is about one-tenth 
of a photon in the cavity, and two atoms interacting with 
the mode. Negative peaks appear at a frequency less than 
that of the coherent coupling, g, between the atoms and 
the field. Figure ^ shows the result of directly solving the 
master equation with the quantum regression theorem. 

The connection between the spectrum of squeezing and 
the time domain fluctuations shows that the field under- 
goes an exponential damping which is not present for 
the weaker driving fields (see Fig. ||ii). The careful reader 
may think there is a slight problem in saying that the con- 
nection between the spectrum of squeezing and the third 
order correlation function still holds in these stronger 
fields. We will show, from the results of a quantum tra- 
jectory simulation in section Vl.b, that the qualitative 
behaviour of the fluctuations remain the same. Therefore 
the third order correlation function provides a legitimate 
method for studying the spectrum of squeezing beyond 
the weak field approximation. 

We can get a preliminary understanding of the strong 
field spectrum of squeezing from the OB work of Reid and 
Walls I^J]. They showed that in the optically bistable 
regime of a strongly driven, on resonance cavity, the 
spectrum of squeezing is a Lorentzian centered at zero 
frequency. We understand that this peak comes from the 
coupling of the cavity mode to its environment, but until 
the recent experiments of Foster et al. |^ , we have had 
no motivation to understand the origin of this peak in 
the time domain. 
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FIG. 8. i. Spectrum of Squeezing and ii. its inverse fourier 
transform for two atom system under a stronger driving field. 
This figure shows the influence of the peak at zero frequency 
on the fluctuations in the time domain. X — 18.1 and 
(.g,K,7)/(27r) = (38.0/72,8.7,3.0) MHz. 

VI. DISCUSSION 

We have analyzed some of the many trajectories which 
go into making up the third order correlation functions 
shown in the last section by setting r « 1 in Eq. ([2^). 
We first begin by reviewing some results in the weak field 
regime and then we consider the more interesting regime 
of stronger driving intensities. 



A. Weak Field Trajectories 

Figure ^ shows examples of the conditional field evolu- 
tion calculated quantum trajectory simulations. Figure 

shows the evolution of the field following the detec- 
tion of a photon escaping through the cavity mode and 
Fig. ||ii shows the field evolution following the sponta- 
neous emission of a photon out the side of the cavity. We 
used the following realistic parameters from the work of 
Foster et al |,|: (5, k, 7)/27r=(38.0, 8.7, 3.0) MHz. This 
means that the system is in the strong coupling regime 
for one and two atoms. 

The basic dynamical mechanism which describes these 
two results can be traced back to the change of the state 
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following a collapse operation of the type found in 
Eq. We now consider the reduction of the equilib- 
rium state of the cavity QED system on the occasion 
of a triggering photon detection (see Fig. p. Defining 
Ag = (aexp(— z0)-|-a^exp(j0))/2, where a is the annihila- 
tion operator for the cavity field and 9 is the BHD phase, 
we consider the quadrature amplitude, Aqo , in phase with 
the steady state of the field A = (a). For weak excita- 
tion, and assuming fixed atomic positions {r}}, to second 
order in A the equilibrium state is the pure state pTUT^ 



IV-Ss) - [|0)+A|l) + (AVA/2)a/3|2) 
+ [</.|0)+A</)/3|l) + ---]|i?) + - 



(30) 



where \G) is the one or two- atom ground state and \E) 
is for one atom in the excited state with all others in the 
ground state. We assume that all the atoms are coupled 
to the cavity mode with the same strength, g. The ex- 
plicit form for 0, a and /3 follow from solving the master 
equation in the steady state They are 



2^Ng 
7 



A 



a 



/3 = 



1 + 2C 



(31) 



(32) 



(33) 



1 + 2C - 2C[ 

where C = NCi and C[ = -^j^^y 

After detecting the escaping photon, the conditional 
state is initially the reduced state a|'0)/A, which then re- 
laxes back to equilibrium. The reduction and regression 
is traced by 0,|l|l 



where 



{|0)+A[l + C/(r)]|l)+---}|G) 



4c;c 



1 + 2C - 2C[ 



(34) 



(35) 



/(r) ==exp(-(K-H7/2)T/2){cosl7T-$sinr2r}, (36) 



2k -I- 7 
AD. 



(37) 



(38) 



Thus, the quadrature amplitude expectation makes the 
transient excursion (Aoo) A[l + C/(''")] away from its 
equilibrium value (Aqo) = A. 



In contrast, an undetected spontaneous emission pro- 
duces the reduced state a~\'ip)/\, which sets up a com- 
pletely different evolution as shown in Fig. ^i. The re- 
duction and regression in this case is again traced by 
Eq. ( ^ ) with the following replacements: 

2C[ 



c = 



2k 



4r2 



1 + 2C- 2C[ 

jn{(3 - 1) 



(39) 



(40) 



N- 



with q — yl 

These two distinct behaviors correspond fairly loosely 
to the regression to equilibrium observed in the step ex- 
citation in the field. Fig. |£i, and a step excitation in the 
atomic polarization. Fig. |ii. Note the phase shift be- 
tween the two responses. The steady state wavefunction 
determines the size of the steps. 



(^0° (cavity)) 
° (steadystate) / 

° (spontaneous) / 



a(3 



(41) 



(42) 



° (steadystate) / 

We see from Eq. (^3|) that /3 is a number greater than 
zero. The field will always remain positive when an atom 
spontaneously emits. As long as g > + 2 )/-^' then 

through a (Eq. 32), the field will change sign when a 
photon escapes through the cavity mode. 
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FIG. 9. Regression of a cavity QED system back to steady 
state; N=l, low intensity i. after the detection of a photon 
escaping out of the cavity mode ii. after the escape of a photon 
through spontaneous emission. The inset shows the sequence 
of events in terms of the cavity QED system and the detector. 
The parameters used are the same as in Fig. ^. 
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The parameters used to create the plots in Fig. ^ are 
the same as those used in Fig. |[ Except for the residual 
shot noise in Fig. these two plots are in agreement 
with each other for t > 0. 



B. High Field Trajectories 

The weak field calculations of section Vl.a make it clear 
that in the strong coupling regime a cavity emission will 
always produce a negative shift in the field. In order for 
us to observe fields of the type in Fig. ^i and spectrum 
of squeezing of the type in Figs. ^, || and ||i we must find 
a mechanism which causes the expectation of the cavity 
field to jump positive after a cavity emission. 

It follows from Eq. ( |30| ) that the ratio of the probability 
for a spontaneous emission to the probability for a cavity 
emission from steady state is 



-fcavity 



(43) 



This implies that in the strong coupling regime it is 
more likely for the system to spontaneously emit out the 
sides of the cavity than directly into the cavity mode. 
Since we only begin a conditioned field measurement af- 
ter a cavity mode emission, then it becomes necessary 
to consider what happens in the likely event of a cavity 
photon following a spontaneous emission. 

Quantum trajectories allow one to follow the dynamics 
of the system under higher excitation. Figure 10 shows 
representative trajectories calculated from Eq. (27) when 
we have two atoms in the cavity and r w 1 . The complete 
trajectory is a series of these excursions from the steady 
state well separated in time, and a returning to the steady 
state as shown. In Fig. |l^ the evolution starts with 
a spontaneous emission (A) out the side of the cavity, 
followed at (B) by a photon escaping through the cavity 
mode that gets registered by the APD detector. Note 
how the field jumps positive and changes curvature with 
the escaping cavity photon. This is exactly the type of 
process which gives rise to the incoherent peak in the 
spectrum of squeezing. 

We offer the following qualitative explanation of these 
types of events. The driving field (e/ft), atom-field cou- 
pling (g), and decay rates (k, 7) are such that the system 
is in a regime where the cavity field is bunched. If we 
have a spontaneous emission event when the system has 
few excitations we return to the steady state value as in 
Fig. If the spontaneous emission event happens while 
in the bunched regime, followed by a cavity emission, 
then there are probably more excitations in the system. 
With one of the atoms removed from the system following 
the spontaneous emission, the probability for this energy 
to be in the cavity mode is increased. If we detect a 
cavity photon soon after the spontaneous emission, then 



we are probably in a regime where the intracavity field 
undergoes a large amplitude fluctuation, and the value 
of the cavity field is higher than the steady state value. 
This causes an upward jump in the expectation of the 
field. This is similar to the upward jump in the con- 
ditioned field of an OPO operated well below threshold 
pO| . In that system the conditioned steady state is small, 
since the system is most likely in the ground state, with 
a small probability of having a pair of photons in the 
cavity. When a photon emerges from the cavity one pho- 
ton remains, and the conditioned average field rises. It 
is also clear that these types of events increase linearly 
with the number of atoms in the cavity, since the ratio 
of spontaneous emission events to cavity loss events is 
2NCi. 

If we drive the same system much harder the time 
evolution of the conditional field shows multiple jumps, 
some from spontaneous emission and some from escapes 
through the mode. The dynamics get very complicated 
and we show an example in Fig. ^^i for illustration, but 
the general trend remains, that the average value of the 
field is much larger than the steady state in such cases. 




12 3 

T/2gT 

FIG. 10. A quantum trajectory simulation which shows the 
time evolution of the field back to steady state after the de- 
tection of a photon, i. A spontaneous emission event followed 
by a cavity emission which starts the averaging process, ii. A 
spontaneous emission event followed by many cavity emission 
events. The inset shows the sequence of events in terms of 
the cavity QED system and the detector. Both figures were 
prepared with the following parameters: N = 2, X = 18.1, 
(p,K,7)/(27r) = (38.0/72,8.7,3.0) Mhz. 

Figure |l^ demonstrates the power of studying the sin- 
gle trajectories and also provides the major result of this 
paper. We see that the field will have a noise background 
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arising from the mean cavity field, on top of other quan- 
tum fluctuations, when one has a spontaneous emission 
event out the side of the cavity. This is usually fol- 
lowed by a series of cavity emissions. This effect depends 
strongly on the number of atoms in the system and also 
on the strength of the driving field. This can not be seen 
from a direct numerical solution of the master equation, 
but now, thanks to the wave-particle correlator unravel- 
ling of the master equation, we begin to see exactly how 
these squeezing spectra are degraded by the incoherent 
process of spontaneous emission. 

The entire trajectory is a collection of events well sep- 
arated in time of the type in Figs. |l^ and |l^i. If we 
average over many random realizations of these different 
events with an initial cavity emission setting the trigger 
at t = 0, then we recover the averaged conditioned field 
evolution. Figure |ll| shows such an average field which 
after symmetrization and upon taking the Fourier trans- 
form leads to a spectrum which resembles the experimen- 
tal observations qualitatively. This also shows a similar 
behaviour as the oscillation from Fig. ^i and therefore a 
legitimate connection can still be made between the cor- 
relation function of Eq. ^ and the spectrum of squeezing 
slightly beyond the weak field regime. 
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FIG. 11. Conditioned field evolution calculated from 5000 
realizations of a photon counting, quantum trajectory, simu- 
lation for N=2 and high intensity. X = 18.1, {g, k,j)/{2tt) — 
(38.0/v^,8.7,3.0) Mhz. 



C. Time Symmetry 

One result that we have only briefly mentioned is the 
fact that Fig. ^ is symmetric with respect to t = while 
the single events presented in Figs. ^ and ^show a clear 
time asymmetry. This symmetry is recovered as a result 
of the BHD back action. The back action comes about 
because of the quantum superposition of both the local 
oscillator and cavity signal field (see Eq. (|2^)). It is not 
correct to think of these two fields as separable, a classical 
notion, but that each time a local oscillator photon is 
detected, it also effects the evolution of the cavity field. 
This implies that the conditional system state at time t is 



correlated with the shot noise that has appeared over the 
recent past. Thus, through the "start" clicks (particle 
aspect), we postselect a subensemble of the shot noise 
which has in effect been filtered through the system's 
dynamical response function (Hamiltonian evolution). 

Quantum trajectories tend to reveal precisely those 
physical attributes they are trying to measure. By adding 
the homodyne detector we attempt to measure the fluc- 
tuations in the field amplitude, and when we do this these 
fluctuations appear in the averaged photocurrent before 
the arrival of a triggering photon. 



D. Future Work 



By stating the problem in the time domain, we open 



the possibilities to apply quantum feedback theory ^6 27 
in order to conditionally stabilize the atom-cavity sys- 
tem evolution against decoherence from the environment. 
This stabilization will not be optimal since the informa- 
tion obtained is limited to the quasimode leaking out of 
the cavity and any spontaneous emission event will be 
missed in the feedback loop. This calls for operating the 
cavity QED system in the low intensity regime to min- 
imize spontaneous emission. We are currently investi- 
gating the theoretical and experimental requirements to 
develop such a program. 



VII. CONCLUSIONS 

The conditional evolution of the field of a cavity QED 
system in the strong coupling regime shows non-classical 
behavior. This is purely quantal in nature and reflects 
the intrinsic dynamics triggered by the collapse of the 
wavefunction after the detection (conditioning) of a pho- 
ton. 

Looking at the time evolution of the conditional state 
as it regresses back to steady state we find that its be- 
havior is far from random. When uninterrupted by other 
events, it follows a dynamic path that shows the coupling 
between the atom and the cavity and the significant ex- 
cursion away from equilibrium brought back by a single 
quantum measurement. The spectrum of squeezing reg- 
isters exactly the same behavior, but without the time 
and phase information critical to understanding the dy- 
namics of the measurement process. 

As we implement the measuring device used in the lab- 
oratory, we require an averaging process to observe the 
conditional evolution of the electromagnetic field. That 
evolution is non-classical. The phase is exactly tt from 
the driving field and it ocurrs when the intensity has a 
fluctuation up. This is contrary to all classical intuition 
that would say that then the field should have also a 
maximum, but in fact the field shows a minimum. 
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The presence of the homodyne measurement intro- 
duces back action into the system that then provides the 
foundation for the time symmetric wave-particle corre- 
lation function when third order correlations can be ne- 
glected. 

The presence of other avenues for the interruption of 

the evolution back to steady state, such as spontaneous 
emission, modify the correlation function, decreasing its 
non-classicality, while at the same time modify the spec- 
trum of squeezing around the origin with a peak that 
shows just the empty cavity width for the one atom case, 
or the spontantous emission width once the system is 
allowed to have more than one atom. 

The wave-particle correlation function opens the door 
to further the study of quantum optical systems in the 
time domain that are necessary to understand better the 
possible realizations of quantum feedback systems. 
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